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Abstract 



We study the rational solutions of the discrete version of Painleve's 
fourth equation (d-Prv)- The solutions are generated by applying 
Schlesinger transformations on the seed solutions —2z and —1/z. After 
studying the structure of these solutions we are able to write them in 
a determinantal form that includes an interesting parameter shift that 
[ — ■ vanishes in the continuous limit. 



1 Introduction 

One important question in the study of discrete versions of continuous dif- 
ferential equations concerns the existence of corresponding special solutions. 
For continuous Painieve equations rational and special function solutions are 
^ known fl], g], and in many cases even a rigorous classification has been done 

0. If one proposes a discrete version of a Painieve equation it is not enough 
that in some continuous limit the continuous Painieve equation is obtained, 
but in addition the proposed equation should share some further properties 
of the original equation. One of these properties should be the equivalent 
of the Painieve property, called "singularity confinement" This has al- 
ready been used to propose discrete forms of the Painieve equations 0. Other 
structures of the continuous Painieve equations that have been shown to exist 
for the discrete ones include their relationships by coalescence limits and 
the existence of Hirota forms for these equations. What is still largely an 
open question is the fate of the special solutions (rational, algebraic, special 
function) known for the continuous case. 
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Here we discuss the determinantal structure of the rational solutions to 
the discrete fourth Painleve equation, (d-Prv), given by 



(x n + a + (3){x n + a - (3)(x n - a + f3)(x r , 



a 



0) 



(x n + z n + i){x n + z n - 7) 

where = <5n + (. [Note that this equation is invariant under the change 
sign of any of parameters a, j3, 7.] 



One reason for calling (|) d-Prv is that if we put 



—a + 



-6/2 



r 2 , /? 2 

Piv: 



-6/2 + 5 2 , 7 = 5 1 and then take the limit 5 — > we get 
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The comparison of ([I]) and (fj) reveals the first interesting difference between 
the discrete and continuous versions: the discrete one has more parameters. 

But why should one insist on writing the solution in determinantal form? 
It is well known that most integrable systems have multisoliton (and rational) 
solutions in determinantal form, and the same holds for many special solu- 
tions of continuous Painleve equations. This determinantal structure actually 
carries fundamental information and reveals the "basic" object hidden in the 
solutions. For example it is well known that Pn, 



d 2 u 
dz 2 



2u 4 - Azu + 4a, 



where a is a parameter, admits rational solutions for a = 
been shown that these solutions can be are expressed as 



u 



d , ttv + i 
— log 
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N + 1. It has also 
1 

(4) 



where 77v's are polynomials in z (Vorob'ev-Yablonski polynomials). However, 
it has been shown only recently that can be expressed in determinantal 
form IH, 



tn 



q N (z,t) 
q N - 2 (z,t) 



q N+1 (z,t) 



q 2 N-i(z,t) 
q 2N - 3 {z,t) 

qi(z,t) 



Q-N+2{z,t) q-N+z{z,t) 
where g^'s are the so called Devisme polynomials defined by 
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(5) 



fc=0 



exp ( zX + t\ 2 + -A 3 



and qk(z,t) = 0, for k < . (6) 



Thus, these Devisme polynomials may be considered as basic polynomials in 
the rational solutions of Pn. Similar determinantal structure of the rational 
solutions is observed for the discrete case: in fact, it has been shown that 
in the case of d-Pn, the rational solutions can be expressed by determinants 



whose entries are given by Laguerre polynomials [10, 11 



In the case of Prv, it is well known JTJ that the continuous Prv has three 
rational solution hierarchies, whose "seed" solutions are 

1 2 

y(z) = -2z, — , or - -z, (7) 
z 3 

and determinantal forms of the first two hierarchies have been discussed in 
general terms, e.g., in 0. The detailed results for the first two hierarchies 
are as follows [0 . Let be an N x N determinant of Hankel type given by 

r v N = det |fr w+i+i _ 2 | 1 < ijJ .<j V , t% = 1 , (8) 

where H n , n = 0, 1, 2 • • • are Hermite polynomials in x characterized by the 
recursion relations 

H n+1 = xH n - nH n _ u -^-# n = nH n _ x , H = 1, H l = x. (9) 

Then 

—^tM^)' z = 7Z x ' (10) 

are rational solutions of Prv© for parameter values 

(a, 6) = (-(V + 2iV+ l),-2z/ 2 ), v,NeZ, v>l, N > 0. (11) 
Moreover, 

-^QMW-*)' ^tt- (12) 

yield rational solutions of Prv(@) with 

(a, 6) = (N- u,-2(u + N + l) 2 ), v,NeZ, v > 0, N>0. (13) 

The solutions given in ( |ToD and (|12D correspond to the "— 1/z" and "—2,2" 
hierarchies, respectively^]. 

Here our object is to find the discrete versions of these results for the 
—2z and —1/z hierarchies. (We hope to return the more complicated —\z 
case elsewhere.) We cannot simply discretize the continuous results, because 



1 Equations ([L0|) and (12) actually give half of the solutions in —1/z and —2z hierarchies 
obtained by Murata ||. The other half is given in terms of polynomials similar to the 
Hermite polynomials, but with some different signs in the recursion relations above. 
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the discrete equation has more parameters, and the way the new parameters 
modify the continuous results is indeed one of the interesting questions. In 
approaching this problem we will not use any detailed properties of the con- 
tinuous case, we just assume that that such rational solutions should arise 
from determinants of polynomials, which always happens in the continuous 
case. Starting with the discrete version of the usual seed solutions, we first 
construct a set of solutions in both hierarchies using Backlund-Schlesinger 
transformations (discussed in the next section) and then search for the deter- 
minantal structure by studying the properties of these rational solutions, the 
main clue being the factorization of the denominator. 



2 Backlund-Schlesinger transformation 

In order to generate rational solutions we use the Backlund-Schlesinger trans- 



formations on the seed solutions. These transformations were given in |L3 
We write them as follows: 

Backlund transformation: Let us assume that x(n) solves d-Prv with 
parameter values (a,/3, 7) (in which case we often write x = x(n; a, (3, 7) or 
x = {fin)] a, (3, 7}). Using these let us define 

xx + x(z + a) + x(z — a) + b 2 — a 2 ) 



BT a:b (x) :-- 



x + x 



(14) 



where z = (n + \)5 + £, and x = x{n + 1). Note the shifts in n. Then we get 
new solutions from the elementary Backlund transformations BT{ as follows: 



x(n; a + 5 5, 7, 0) = BTix{n) = BT a ^(x{n; a, (3, 7), 

x(n; a - |5,7,/3) = BT 2 x(n) = BT_ a ^(x(n] a, (3, 7), , . 

x(n; (3 + 5 5, 7, a) = BT 3 x(n) = BTp >a (x(n; a,/3,j), 

x(n; (3 — |5, 7, a) = BT A x{n) = BT_p >a (x(n] a, f3, 7). 

These transformations jump too much in the parameter space and therefore it 
is useful to define Schlesinger transformations that generate new solutions 
by changing only one of the parameters 

x(n] a + 5, (3, 7) = Si(x) = BT\BTi x(n; a, (3, 7), 

x(n; a - 5,(3, 7) = S 2 (x) = BT 2 BT 2 x(n; a, (3,^), , , 

x(n; a, (3 + 5, 7) = S 3 (x) = BT 2 BT^BTiBT^ x(n; a,(3,j), 

x(n; a, (3 — 8, 7) = Ss(x) = BT 2 BT%BT 2 BT& x(n; a, (3, 7). 

It should be noted that sometimes there are barriers over which Schlesinger 
transformations cannot cross. For example ifa = 7+(m + |)5 and (3 = 
7 — (n + |)5 (which is relevant for some rational hierarchies) then the barrier 
n = cannot be crossed by Schlesinger transformations because BT 3 yields 
the solution, similarly for m = with BT 2 . 
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3 Some solutions 

Using the above Schlesinger transformations we can construct other rational 
solutions from the seed solutions. The seed for the — 2z hierarchy is given by 

£oo = {-2z,7 + i5,7-^,7}. ( 17 ) 

(here and in the following z = n5 + (.) Since equation (|1]) is invariant under 
the sign changes in a, and 7 we may always assume the above signs in front 
of 7. With Schlesinger transformations we can reach the parameter values 



a 



7 + (M+i)5, /3 = 7 -(jV+i)5, 



where M and N are nonnegative integers, and we may also assume that 
M > N because of the a <-> (3 symmetry. This then defines the extent of the 
hierarchy. For particular solutions we give N, M as subscripts. Some further 
solutions obtained this way are (see also fT4|| ) 

= < y ' 7 + i 5 ' 7- i 5 ' 7 f' 



1 4z 3 + 2^(-3 7 -55) 1 . 1 

X20 = < -2^ + 5(7 + 25) ' 7 + ^7-^7^, 



/ 4z 4 + 4r^(-3 7 - 75) + 3<5 2 ( 7 2 + 5 7 £ + 65 2 ) 7 1 1 
X3 ° ~ j -2^ + z5(3 7 + 85) >7+5 d »7 2^7j, 

/ 8z 5 - 16z 3 ^ + 2z^ 2 (3 7 2 + ^) 3 _3 1 
X31 ~ j -4z 4 + 4z 2 5 2 + 5 2 ( 7 2 -5 2 ) ' 7+ 2 d '7 ^7 j • 

The — 1 /^-hierarchy is not connected to the — 2z hierarchy by a Schlesinger 
transformation, but only by a Backlund transformation followed by a redefi- 
nition of 7. In any case the seed is 

xio = { "^ +<y) , 7+l*,7 + i*}, (19) 

and in general the parameter values in this hierarchy are 

a = 7 + (M + i)5, /? = 7 + (jV + i)5, (20) 
with < iV < M. Some further solutions are given by 
j 2z5(2-f + 35) 



X20 



, 7 + §£, 7 + h 5 , 7 ^ , 



-2^ 2 + 5(7 + 25) ' 

<5(7 + 25)(-2^ 2 + <5(7 + 5)) . , s _ , 



X21 = j 2^5^ ,7 + ^7+2^7 



tffo 3 7 + 1^ 

L I Z U [3][2] [4][6] Z 1 ^ 2 U 

J2]_ _[5]_ J8L 1 ,3r 

U [1][2] [2][4] [3][6] 1 >" 

_ LU [2] [3] 

0][1] [0][2] [0][3] [0][4] 



J1L_ J2i_ J1L n ^ + 1^ 

u n i mim n a n u u l * ■?. u 



M= 12 3 4 

a= l+¥ 7 + §5 7 + \5 7 + |5 



N /3 



n _t°L _IU_ _LM_ n 

U [0][1] [0][2] [0][3] [0][4] U 7 2 

[3] [6] [9] [12] 1 7 _3^ 



[2][2] [3][4] [4][6] [5][8] x ' 2 

-2z 2 

ZZ [6][6] [9][8] / f 2° 



[231 



[12][12] 



3 -y-ls 



Figure 1: Observed degrees in n (in square brackets) of various solutions. The 
denominator always factorizes. For the — 2z hierarchy we have subtracted — 2z 
from each solution. 



x 30 = j 3*(-7-2*)(-2^(7 + 3^)) 

^ -2z 3 + zd(3 1 + 85) ,7+ 2 o,7+ 2 o,7 



j (_ 7 _ 35)5(4z 4 - AzH* + 3£ 2 ( 7 2 + 2 7 £ + 5")) \ 
X32 \ 4z 5 + 4z 3 5(2 7 + 5) +^ 2 (3 7 2 + 4 7 5 + 5 2 ) ' 7 2 ' 7 2 ' 7 J ' 

It is obvious from these examples (and from the equation) that there is an 
overall scaling invariance, and because of this we will in the following simplify 
expressions by scaling out S by writing 7 = cS, z = nS. 



4 The elementary polynomials 

A common property for rational and other special function solutions of the 
continuous case is that the denominator factorizes into two determinants. 
Thus after constructing a set of solutions (all computations were done using 
REDUCE Hl5| ) we studied the factorization of their denominators, the results 
are given in Figure 1. From it we can see that for the — 2z hierarchy the 
denominators of xmn factorize with factors of degree N(M+1) and (N+1)M 
in n, and for the —1/z hierarchy the degrees are N(M — N) and (N + 1)(M — 
N). 
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Next observe that when N = one of the factors in the denominator 
is constant and the other grows linearly with M. (For this value N = 
the denominators are the same for the — 2z and —1/z hierarchies.) The 
denominator is interpreted as the product of an x and lxl matrix, and 
this suggests that the denominators can provide us the basic polynomials that 
are to be used as matrix elements. The first few polynomials obtained this 
way are (when normalized to be monic): 



Po 


J 1 , 


c) 


■■= 1, 




Pi 


>, 


c) 


:= n, 




P2 


>, 


c) 


:= n 2 — \c — 1, 




P3 




c) 


:= n 3 + n(— |c — 


4), 


Pi 


% 


c) 


:= n 4 + n 2 {— 3c - 


-10) + |c 2 + fc + 9, 


P5 


% 


o) 


:= n 5 + 5n 3 (-c- 


-4)+n(fc 2 + ifc + 64), 


P6 


% 


c) 


:= n 6 + n 4 (-fc 


-35)+n 2 (fc 2 + ^fc + 259) 








_15 3 _ 225 

8 8 


c 2 _ 555 c _ 225; 


P7 


% 


o) 


:= n 7 + n 5 (-fc 


-56)+n 3 (^c 2 + ±±pc + 784) 








+n(- 1 f-c 3 


_ 1|5 C 2_2|9 C _ 2304) 



One indication that we are on the right track is obtained when we observe 
that the p^{n 1 c) satisfy recursion relations: 

p N+1 (n,c) = np N (n,c) - f (c + N + l)pjv-i(n, c + 1), ^ 
PN{n + |, c) - p N (n - |, c) = Np N _!(n, c + \). ~ 

As was noted in Sec. 1 the matrix elements of the corresponding continuous 
case are given in terms of Hermite polynomials. Now we note that in the 
continuous limit (which means n = x /(y/26), c = 1/6 2 , 5 -> 0) we get for 

H N (x) := lim (V26) N p N (x/(V2S), 1/5 2 ), (22) 
the recursions relations of Hermite polynomials ([|). 



5 Matrix form of the denominators 

The above indicates what the polynomial matrix entries should be, and de- 
termining the structure of these matrices in the denominator is the next 
problem. 

For N = the denominator was interpreted as the product of an x 
and lxl matrix. For iV = 1 we can see from the table that one of the 
factors grows by 1 as M increases by 1 and the other factor increases by 2. 
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These factors were interpreted as 1 x 1 and 2x2 matrices, respectively. The 
lxl part for N = 1 was then found to be proportional to the shifted basic 
polynomials given above: for a given M this factor is pm+i(c — 1) in the — 2z 
hierarchy and £>m+i(c+ 1) i n the — 1/z hierarchy. The necessity of shifts in c 
is an important new ingredient, and something that exists only in the discrete 
version. 

The next problem was to find a suitable matrix structure for the factors of 
the denominator. Our working assumption was that for any M the denomi- 
nator should be a product of a A x A and a (A + 1) x (A + 1) matrix. The 
degrees of these factors led us to try the determinantal structure of Hankel 
type 



L,K 



(lK+l 



Qk+i 

QK+2 



Qk+l-i 
Qk+l 

QK+2L-2 



Qk+l-i Qk+l 

because now if the entries qx are non-monic polynomials of degree K then 
deg{T L>K )=L(K + L-l). 

A more detailed study with some trial and error revealed that q's are not 
simply proportional to the p's obtained before, but that different shifts in c 
are also needed for the matrix elements. Our final result is as follows: 

For the -2z hierarchy (for which a = 7 + (M + \)5, (3 = 7 - (N + \)5 
where < iV < M) the denominator can be expressed as 



den(x M ,N) = t NjM -n+2(-1) t n+1 , M -n(0), 



(23) 



where 



tl,k( s ) '■ = 



Qk(s) 
q K +i(s - 1) 



Qk+i(s) 
q K +2(s - 1) 



qK+L-i(s - L + 1) q K+L (s-L + l) 



and 



q K +L-i(s) 
q K +L(s - 1) 

qx+2L-2(s - L + 1) 



,(24) 



q M (s) := ( c + s) M p M (c + s)/M\. (25) 

For the -1/z hierarchy (a = 7 + (M + 1)5, (3 = 7 + (N + \)5, where 
< < M) we found same matrix form, but with different degrees and 
shifts: 

den(x M ,N) = r NjM -2N+i(N) t n+1iM _ 2 n(N). (26) 
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6 The numerator 



Finding a determinantal form for the numerator was more difficult, because 
it did not factorize. However this was expected, because usually it turns out 
that the numerator is the sum of two products of r-functions, c.f., (p|,|iO|,|T2h . 

Thus we tried to express the numerator as a sum of two products of two r- 
functions (£23|), with possible shifts not only in c, but also in n, corresponding 
to the derivatives in the continuous cases. Furthermore it seemed reasonable 
to assume that the sizes of the matrices were the same as in the denominators, 
and that the other index of the r-function depended linearly on M and N. 
Figuring out the shifts in n required some more trial and error, but eventually 
we arrived at the result that worked: 

For the — 2z hierarchy we got 



x{n- 7 + (M + \)8, 7 - (N + §)<$, 7) oc 

COsll(|.D n ) T7v,M-AT+l( — §) - TN+l,M-N+l( — §) 
Tn,M-N+2( — 1) 7"JV+1,M-Jv(0) 

and the corresponding result for the —1/z hierarchy was 

x(n; 7 + (M + 1)5, 7 + (N + \)8, 7) oc 

sinh(|D n ) t N)M -2N+i{N - |) ■ r Af+ i i A f _ 2 Af(A^ + \) 

T~N,M-2N+l{N) Tn+1,M-2n(N) 



(27) 



(28) 



In expressing the n-shifts we have used the usual Hirota bilinear deriva- 
tive operator D n . These results look a bit more symmetrical if we define 
tn,a+b+n-i{2B - N + 1) := t NjA (B), then we get 



x(n; 7 + (M + \)8, 7 - (N + |)<5, 7) oc 

cosh(|D n ) f N<M _i(-N) ■ f N+1M+ i(-N - 1) 
tnm(- n ~ !) t n +i,m(-N) 

and 

x{n- 7 + (M + |)5, 7 + (N + 1)8, 7) oc 

sinhd^) f N>M _i{N) ■ f N+1>M+ i{N + 1) 



(29) 



T N,M 

(N) 



(30) 
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7 Conclusions 



We have shown here that the two hierarchies of rational solutions for d-Prv 
(generated from —2z and —1/z) can be expressed in terms of determinants, 
and that the matrix elements of these determinants are given by polynomials 
that can be regarded as discrete analogues of Hermite polynomials. It is 
interesting that in these expressions the parameter c, which vanishes in the 
continuous limit, plays an important role. However, the result is not yet 
complete, because the bilinear equation for the r-function is still to be written, 
and the general proof must be given. 

For the —\z hierarchy things are an order of magnitude more difficult. 
There is no linear growth in any direction in the parameter space so there are 
no candidates for matrix elements. We have nevertheless found what the r- 
functions should be, but no determinantal expression for them. Although no 
determinantal expression or basic polynomial is known even in the continuous 
case, we hope that similar structures in solutions of both continuous and 
discrete cases will be found. 
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